Berry phase effect in anomalous thermoelectric transport 
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We develop a theory of Berry phase effect in anomalous transport in ferromagnets driven by sta- 
tistical forces such as the gradient of temperature or chemical potential. Here a charge Hall current 
arises from the Berry phase correction to the orbital magnetization rather than from the anomalous 
velocity which does not exist in the absence of a mechanical force. A finite-temperature formula for 
the orbital magnetization is derived, which enables us to provide an explicit expression for the off- 
diagonal thermoelectric conductivity, to establish the Mott relation between the anomalous Nernst 
and Hall effects, and to reaffirm the Onsager relations between reciprocal thermoelectric conductiv- 
ities. A first-principles evaluation of our expression is carried out for the material CuCr2Se4-a:Bri;, 
obtaining quantitative agreement with a recent experiment. 

PACS numbers: 72.15.Jf,75.47.-m,75.20.-g 



The phenomena of transport fall into two categories: 
those due to a mechanical force, such as the electric field 
on charges, and those driven by a statistical force, such as 
the gradient of temperature or chemical potential. The 
mechanical force exists on the microscopic level and can 
be described by a perturbation to the Hamiltonian for 
the carriers, while the statistical force manifests on the 
macroscopic level and makes sense only through the sta- 
tistical distribution of the carriers. It has been estab- 
lished [l|, I3] that Berry phase of Bloch states has a pro- 
found effect on transport driven by a mechanical force. 
This is through the mechanism that the group velocity 
of a Bloch electron acquires an anomalous term propor- 
tional to the mechanical force, i.e.. 



1 g£n(fc) 
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-E X f2„(fe) , 
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dk h 

where e„(A;) is the band energy, —eE is the mechani- 
cal force due to the external electric field, and r2„(fc) 
is the Berry curvature, the Berry phase per unit area 
in the fc-space. Evaluation of the Hall current from the 
anomalous term reproduces the Karplus-Luttinger for- 
mula [Sj for the anomalous Hall conductivity. Calcula- 
tions based on the Berry phase effect have found much 
success in explaining anomalous Hall effects (AHE) in 
ferromagnets of semiconductors |J| , oxides pj and tran- 
sition metals p . Recent experiments |2l |S| give further 
convincing evidence in support of this theory. 

A natural question is whether and how the Berry phase 
also manifests in transport driven by a statistical force. 
On the one hand, the anomalous velocity vanishes in 
the absence of a mechanical force, eliminating the ob- 
vious cause for a Berry phase effect in this case. On the 
other hand, this conclusion would introduce a number 
of basic contradictions to the standard transport the- 
ory. First, a chemical potential gradient would be dis- 
tinct from the electrical force, violating the basis for the 
Einstein relation for transport. Second, a temperature 



gradient would not induce an intrinsic charge Hall cur- 
rent, violating the Mott relation [see Eq. (|HJ|I below] be- 
tween the AHE and the anomalous Nernst effect (ANE), 
where a transverse current is produced by a temperature 
gradient in ferromagnets. Third, as will be made clear 
below, it would be impossible to establish the Onsager 
relation between cross transport coefficients connecting 
thermoelectric Hall currents and forces. In addition, a 
recent experiment on the ANE in the spinel ferromagnet 
CuCr2Se4_a;Br2, [3 found weak dependence on scatter- 
ing, suggesting that there should indeed be a Berry-phase 
induced intrinsic mechanism. 

In this Letter, we solve the puzzle by showing how the 
Berry phase effect manifests in thermoelectric transport 
driven by a statistical force. It turns out that the local 
current of carriers acquires an extra term from the car- 
rier magnetic moment in the presence of a non-uniform 
distribution which arises from the gradient of tempera- 
ture or chemical potential. However, the complete the- 
ory also relies on a proper deduction of magnetization 
current |lC|, and requires a deeper understanding of the 
orbital magnetization. It was found that there is a Berry- 
phase correction to the magnetization [ill . Il2l | , and here 
we generalize it to the case of finite temperatures which 
is needed for thermoelectric transport. This Berry phase 
correction eventually enters into the transport current 
produced by the statistical force, playing the counter- 
part as the anomalous velocity term due to a mechanical 
force. 

We have thus found perfect harmony between sta- 
tistical and mechanical forces even in the presence of 
Berry phase effect. The basic transport relations of Ein- 
stein, Mott, and Onsager continue to hold, which gives 
strong support for the validity of our theory. Finally, 
we also provide a reality check on the Berry phase ef- 
fect in the ANE by calculating the intrinsic anomalous 
Nernst conductivity axy ll3'| for CuCr2Se4_2:Bra; using 




wave packet 
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FIG. 1: The wave packet description of a charge carrier whose 
center is {vc, kc). A wave packet generally possesses two kinds 
of motion: the center of mass motion and the self-rotation 
around its center. Both of them contribute to the local current 
density as given in Eq. (|5J. 



first-principles method. The obtained doping dependence 
curve agrees well with available experimental data |3|. 
Our calculation also predicts a peak-valley structure be- 
tween the data points, at a place where the anoma- 
lous Hall conductivity has a sudden sign and magnitude 
change. 

Local and transport currents. — In the conventional 
Boltzmann transport theory, one considers a statistical 
distribution g(r, fc) of carriers in the phase space of po- 
sition and crystal momentum. The distribution function 
satisfies the Boltzmann equation with a collision integral 
whose form depends on the details of the collision process. 
The current density is given by J = —ej[dk]g{r,k)r, 
where J[dk] is a shorthand for J dk/{2TT)^, and a sum- 
mation over band index has been omitted for simple no- 
tation. In the absence of a mechanical force, the electron 
velocity is simply r = de{k)/hdk. It is then apparent 
that the anomalous velocity term due to the Berry phase 
drops out of the expression for the current. 

However, the above picture is naive that the carrier 
is treated as a structureless point particle. The quan- 
tum representation of the carrier is in fact a wave packet, 
which has a finite spread in the phase space. The wave 
packet generally rotates about its center position, as il- 
lustrated in Fig. ^ giving rise to an orbital magnetic 
moment m(fc) = ~{e/2){W\{f — Tc) x ■w|VF}, wherejVF) 
is the wave packet and v is the velocity operator [l|,|2|- A 
careful coarse graining analysis 14] shows that the cor- 
rect expression for the local current has an extra term: 



: l[dk]g{r,k)r + V X f [dk] f {r , k)m{k) 



(2) 



where the magnetic moment enters explicitly. In the ex- 
tra term we have replaced g{r, k) with the local equilib- 
rium Fermi-Dirac distribution /(r, k) for a linear-order 
calculation. When the temperature or chemical poten- 
tial varies in space, the extra term will be proportional 
to the gradient of these thermodynamic quantities and is 
therefore non-negligible. 

For transport studies, it is important to discount the 
contribution from the magnetization current, a point 
which has attracted much discussion in the past. It was 
argued that the magnetization current cannot be mea- 



sured by conventional transport experiments (For a re- 
cent most comprehensive work, see Ref. [lOj). Therefore, 
one introduces the concept of transport current, defined 
by 



j = J - V X M{r) 



(3) 



where M{r) is the magnetization density. This is entirely 
analogous to the classic distinction between microscopic 
current and free current [l5l |. 

It is also important to realize that the magnetization 
density is not simply a statistical sum of the carrier mag- 
netic moments. It has been shown recently that there is 
a Berry phase correction to the magnetization |ll[ Il2l | . 
The contribution from the carrier magnetic moments to 
the local current will be subtracted out in the transport 
current, but the Berry phase correction to the magne- 
tization will give rise to an extra term in the transport 
current. Earlier work concentrate on the zero tempera- 
ture magnetization, and we provide an extension to the 
finite temperature case below. Using Eq. © for the mag- 
netization, we find that the transport current is given by 



dk]g{r, k)r 



'^ X 4 [[dk]^fl{k)\og{l + e-'^('-^'>) 
P J n 



(4) 



where /3 = l/ksT, and the Berry curvature is defined 
by Jl(fe) = Vfe X {u\i'Vk\u) with \u) being the periodic 
amplitude of the Bloch wave. 

The above expression gives a complete account of the 
transport current in ferromagnets, and for crystals with 
nonzero Berry curvatures in general. The first term is the 
usual expression for the charge current, which vanishes at 
local equilibrium (assuming the absence of a mechanical 
force), i.e., g{r,k) = f{r,k). Nonequilibrium correction 
to first order in the gradient of temperature or chemi- 
cal potential yields a result strongly depending on the 
relaxation process, and a transverse current can result 
from skew scattering due to spin-orbit coupling [l6j . The 
second term is new, which results from the Berry phase 
correction to the magnetization. It is also first order in 
the statistical force, but is independent of the relaxation 
time, and is therefore an intrinsic property of the system. 

Orbital magnetization at finite temperatures. — The or- 
bital magnetization of Bloch electrons has been an out- 
standing problem in solid state physics. Recently, dif- 
ferent approaches |llL Il2| have been used to derive a 
formula at zero temperature, where Berry phase is found 
to play an important role. In order to study thermoelec- 
tric transport, we need to generalize it to finite temper- 
atures. Our derivation is made easy by usin g th e field- 
dependent density of states introduced in Ref. [ij , where 
it was shown that in the weak-field limit, a quantum-state 
summation J^k ^i^) of some physical quantity 0{k) 



should be converted to a fc-space integral according to 
J[dk]{l + eB-n/n)0{k). 

The equilibrium magnetization density can be ob- 
tained from the grand canonical potential, which, within 
first order in the magnetic field, may be written as 



F 
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/3 
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P J ri 



(5) 



where the electron energy e^[ = e{k) — m{k) ■ B in- 
cludes a correction due to the orbital magnetic mo- 
ment m(k). The magnetization is then the field deriva- 
tive at fixed temperature and chemical potential, M — 
-{dF/dB)^,T, with the result 



M{r) = / [dk]f{r,k)m{k) 



[dk]-n{k)\og{l + e 



-/3(e-M) 



(6) 



)• 



For generality, we have included a position dependence 
to cover the situation of local equilibrium with a position 
dependent temperature and chemical potential. 

We have thus derived a general expression for the equi- 
librium orbital magnetization density, valid at zero mag- 
netic field but at arbitrary temperatures. The first term 
is just a statistical sum of the orbital magnetic moments 
of the carriers originating from self-rotation of the carrier 
wavepackets. It has been derived in Ref. [J, |2] with the 
expression Tn(fc) = —i{e/2h){'Vku\x[H{k)—e{k)]\'Vku), 
where H{k) is the crystal Hamiltonian acting on \u). It 
has the same symmetry properties as the Berry curva- 
ture. The second term of Eq. © is the Berry phase 
correction to the orbital magnetization. It is of topolog- 
ical nature, arising from a bulk consideration on the one 
hand as in the above derivation, and being connected to 
a boundary current circulation on the other 17]. Inter- 
estingly, it is this second term that eventually enters the 
transport current. 

Anomalous thermoelectric transport. — With the aid of 
Eq. Q it is straightforward to calculate various ther- 
moelectric response to statistical forces. For exam- 
ple, a chemical potential gradient will produce, through 
the second term, a Hall current given by — V/i x 
{e/h) J[dk]f{k)fl{k). This is the same as the Berry- 
phase induced anomalous Hall current in response to an 
electric field if one substitutes V/i/e for the field. It is 
gratifying to see that the Einstein relation continues to 
hold in the presence of the Berry phase effect. 

In the presence of a temperature gradient, an intrinsic 
Hall current also results from the second term of Eq. Q), 



One can then extract an anomalous Nernst conductivity 
axy defined by jx = axy{—'^yT). On a different route, 
we can also obtain the same result by invoking a fictitious 
gravitational field jla| , establishing the Einstein relation 
between this mechanical force and the temperature gra- 
dient. 

Interestingly, by integration by parts, axy can be writ- 
ten into the following more suggestive form 



'^xy 



ds^cFxy[e)- 



T 



(8) 



where cTxy{s) is the intrinsic anomalous Hall conductivity 
at zero temperature with Fermi energy e, given by 



(Jxy{e) 



f[dk]e(e-ek)n,{k). 



(9) 



At low temperatures, the above relation reduces to 



klT 



c^L(£-f) 



(10) 



Such relations between the electrical and thermoelec- 
tric conductivities are known as Mott relations. They 
were proved for non-magnetic materials without or with 
a magnetic field _1£, J£] . Our result extends the validity 
of this relation to ferromagnets and other systems with 
a Berry curvature, and justifies the usage of Eq. pOfl in 
Ref. 01. 

The reciprocal of the ANE is the generation of a trans- 
verse heat current by an electric field. Onsager relation 
dictates that the Berry phase should also affect the lat- 
ter. To show this explicitly, we consider the energy cur- 
rent carried by a wave packet {W\{Hr + rH)/2\W) = 
er — E X Tn{k), where the second term is from the field 
correction to the local Hamiltonian. Assuming a uniform 
temperature and chemical potential |21|, we obtain the 
local energy current to first order in the electric field: 



dk]g{k)er- E X [dk] f {k)m{k) , (11) 



J' 



where the electron velocity r is given by Eq. (Q. How- 
ever, the energy current also has a magnetization part 
from an "energy" magnetization 10] . In the present 
case, it is given by —E x M, which is nothing but the 
material-dependent part of the Poynting vector E x H 
describing the energy flow (with H ~ B/iiq — M) |l 5| . 
Since this energy flow exists in an equilibrium state, it 
does not correspond to a transport current thus must be 
subtracted from J^ to yield the transport energy cur- 
rent j^ = J^ + E X M. Based on our expression (|SJ| for 
the magnetization density, we finally find the Berry phase 
correction to the heat current (defined hy j'^ = j^ — fij): 
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J[dkM{e^fi)f 
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(7) 



Ex-J[dk]n[{s-fi)f 
+ fcBriog(l + e-'3(^-'')) 



(12) 



< 

O 

H 
>■ 
X 




0.4 0.6 

doping X 



FIG. 2: The intrinsic anomalous Nernst conductivity Oxy (di- 
vide by the temperature T) of CuCr2Se4_a;Br2, as a function 
of the Br content x. The calculated curve is compared with 
experimental results • extracted from Ref. Q. 



while the usual expression for the heat current is 
J[dk] g{k){e — fi)v, where v is the usual group velocity 
determined by the band energy. In this case, the Berry 
phase correction comes from both the anomalous velocity 
and the orbital magnetization. Comparison with Eq. ((TJ 
shows that the Onsager relation is indeed satisfied, pro- 
viding a strong evidence for the validity of our theory. 

Comparison with experiment. — The intrinsic anoma- 
lous Nernst conductivity axy only depends on the band 
structure and Berry curvature, so it can be evaluated for 
crystals based on first principles methods. Here we re- 
port our result for CuCr2Se4_a;Br2; and compare with the 
experiment 9] . The band structure and Berry curvature 
are calculated following the procedures in Ref. ^j , using 
the generalized gradient approximation for the exchange- 
correlation potential. Such calculations are very exten- 
sive, and, to reduce the work load, we assume that dop- 
ing affects the Fermi energy but not the band structure, 
which is justified for the present compounds 22]. 

The calculated axy is plotted in Fig. [3 as a function 
of doping X together with the experimental data from 
Ref. JQj]. The comparison is quantitatively good, except 
for the data point at a; = 0.25. This is however a rather 
special point, because it was reported |3| that, for un- 
known reasons, a^y is not really proportional to T for 
X = 0.25. At low temperatures, a proportional rela- 
tion is expected from the Mott relation, which is followed 
strictly by all the data points at other doping densities. 

We also note that while our theory predicts a pro- 
nounced peak-valley structure around x — 0.3, the avail- 
able experimental data at present is too sparse to confirm 
or disprove it. The oscillatory behavior results from the 
complicated band structure of this material, and occurs 
when the Fermi energy (which depends on doping) goes 
through a region of spin-orbit induced energy gap. De- 



tailed explanation based on the numerical calculations 
will be presented elsewhere |22l |. An indirect experimen- 
tal evidence for this peak is that it occurs at a place where 
the anomalous Hall conductivity has a sudden change of 
sign and magnitude around x = 0.3 according to Ref. JT"!. 
Such a correlation is expected from the Mott relation l(Tn|) 
and the fact that the Fermi energy ch ang es approxi- 
mately linearly with the doping density |22]. Neverthe- 
less, more direct experimental results are clearly needed 
for a careful comparison with our theory. 
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